Auto-phylo v2 and auto-phylo-pipeliner: building advanced, flexible, and reusable pipelines for phylogenetic inferences, estimation of variability levels and identification of positively selected amino acid sites.

The vast amount of genome sequence data that is available, and that is predicted to drastically increase in the near future, can only be efficiently dealt with by building automated pipelines. Indeed, the Earth Biogenome Project will produce high-quality reference genome sequences for all 1.8 million named living eukaryote species, providing unprecedented insight into the evolution of genes and gene families, and thus on biological issues. Here, new modules for gene annotation, further BLAST search algorithms, further multiple sequence alignment methods, the adding of reference sequences, further tree rooting methods, the estimation of rates of synonymous and nonsynonymous substitutions, and the identification of positively selected amino acid sites, have been added to auto-phylo (version 2), a recently developed software to address biological problems using phylogenetic inferences. Additionally, we present auto-phylo-pipeliner, a graphical user interface application that further facilitates the creation and running of auto-phylo pipelines. Inferences on S-RNase specificity, are critical for both cross-based breeding and for the establishment of pollination requirements. Therefore, as a test case, we develop an auto-phylo pipeline to identify amino acid sites under positive selection, that are, in principle, those determining S-RNase specificity, starting from both non-annotated Prunus genomes and sequences available in public databases.

Rare serotype c Haemophilus influenzae invasive isolate: characterization of the first case in Portugal.

We report for the first time in Portugal a serotype c Haemophilus influenzae isolated from an adult, with HIV-1 infection. Whole-genome sequencing characterized the isolate as clonal complex ST-7, albeit with a novel MLST (ST2754) due to a unique atpG profile. Integration of this genome with other available H. influenzae serotype c genomes from PubMLST revealed its overall genetic distinctiveness, with the closest related isolate being identified in France in 2020. This surveillance study, involving collaboration among hospitals and reference laboratory, successfully contributed to the identification and characterization of this rare serotype.

Epidemiological Data and Antimicrobial Resistance of Campylobacter spp. in Portugal from 13 Years of Surveillance.

This study extensively analyzed campylobacteriosis surveillance in Portugal from 2009 to 2021, aiming to investigate demographic shifts, seasonal variations, and antimicrobial resistance (AMR) within Campylobacter isolates. Surveillance network and sentinel laboratory-based system data revealed a substantial under-notification of campylobacteriosis cases, suggesting an underestimated disease burden. Notification rates exhibited a paradigm shift, with a notable prevalence among the pediatric population, particularly in children aged 1-4 years, diverging from European reports. Additionally, an emerging trend of Campylobacter infections in younger adults (15-44 years) was observed. The study unveiled a unique seasonal distribution of cases, defying typical summer peaks seen elsewhere. AMR analysis revealed high resistance to ciprofloxacin and tetracycline, in both C. jejuni (93.7% and 79.2%, respectively) and C. coli (96.5% and 93.2%, respectively), stable throughout the studied period (2013-2021). C. coli exhibited significantly higher resistance to erythromycin, gentamicin, ampicillin and ertapenem compared to C. jejuni (p < 0.001). Multilocus Sequence Typing (MLST) data demonstrated the distribution of resistance markers across diverse sequence types, challenging the notion of a clonal origin for multidrug-resistant isolates. In conclusion, the study highlights the need for enhanced surveillance and raises concerns about alarming AMR levels, recommending the implementation of whole-genome sequencing (WGS)-based surveillance for a deeper comprehension of disease patterns and an evolving AMR landscape.

Genome-Scale Characterization of Mycobacterium abscessus Complex Isolates from Portugal.

The Mycobacterium abscessus complex (MABC) is an emerging, difficult to treat, multidrug-resistant nontuberculous mycobacteria responsible for a wide spectrum of infections and associated with an increasing number of cases worldwide. Dominant circulating clones (DCCs) of MABC have been genetically identified as groups of strains associated with higher prevalence, higher levels of antimicrobial resistance, and worse clinical outcomes. To date, little is known about the genomic characteristics of MABC species circulating in Portugal. Here, we examined the genetic diversity and antimicrobial resistance profiles of 30 MABC strains isolated between 2014 and 2022 in Portugal. The genetic diversity of circulating MABC strains was assessed through a gene-by-gene approach (wgMLST), allowing their subspecies differentiation and the classification of isolates into DCCs. Antimicrobial resistance profiles were defined using phenotypic, molecular, and genomic approaches. The majority of isolates were resistant to at least two antimicrobials, although a poor correlation between phenotype and genotype data was observed. Portuguese genomes were highly diverse, and data suggest the existence of MABC lineages with potential international circulation or cross-border transmission. This study highlights the genetic diversity and antimicrobial resistance profile of circulating MABC isolates in Portugal while representing the first step towards the implementation of a genomic-based surveillance system for MABC at the Portuguese NIH.

Multisystemic RFC1-Related Disorder: Expanding the Phenotype Beyond Cerebellar Ataxia, Neuropathy, and Vestibular Areflexia Syndrome.

Background and Objectives: The RFC1 spectrum has become considerably expanded as multisystemic features beyond the triad of cerebellar ataxia, neuropathy, and vestibular areflexia syndrome (CANVAS) have started to be unveiled, although many still require clinical replication. Here, we aimed to clinically characterize a cohort of RFC1-positive patients by addressing both classic and multisystemic features. In a second part of this study, we prospectively assessed small nerve fibers (SNF) and autonomic function in a subset of these RFC1-related patients. Methods: We retrospectively enrolled 67 RFC1-positive patients from multiple neurologic centers in Portugal. All patients underwent full neurologic and vestibular evaluation, as well as neuroimaging and neurophysiologic studies. For SNF and autonomic testing (n = 15), we performed skin biopsies, quantitative sensory testing, sudoscan, sympathetic skin response, heart rate deep breathing, and tilt test. Results: Multisystemic features beyond CANVAS were present in 82% of the patients, mainly chronic cough (66%) and dysautonomia (43%). Other features included motor neuron (MN) affection and motor neuropathy (18%), hyperkinetic movement disorders (16%), sleep apnea (6%), REM and non-REM sleep disorders (5%), and cranial neuropathy (5%). Ten patients reported an inverse association between cough and ataxia severity. A very severe epidermal denervation was found in skin biopsies of all patients. Autonomic dysfunction comprised cardiovascular (67%), cardiovagal (54%), and/or sudomotor (50%) systems. Discussion: The presence of MN involvement, motor neuropathy, small fiber neuropathy, or extrapyramidal signs should not preclude RFC1 testing in cases of sensory neuronopathy. Indeed, the RFC1 spectrum can overlap not only with multiple system atrophy but also with hereditary motor and sensory neuropathy, hereditary sensory and autonomic neuropathy, and feeding dystonia phenotypes. Some clinical-paraclinical dissociations can pose diagnostic challenges, namely large and small fiber neuropathy and sudomotor dysfunction which are usually subclinical.

Viral genetic clustering and transmission dynamics of the 2022 mpox outbreak in Portugal.

Pathogen genome sequencing during epidemics enhances our ability to identify and understand suspected clusters and investigate their relationships. Here, we combine genomic and epidemiological data of the 2022 mpox outbreak to better understand early viral spread, diversification and transmission dynamics. By sequencing 52% of the confirmed cases in Portugal, we identified the mpox virus sublineages with the highest impact on case numbers and fitted them into a global context, finding evidence that several international sublineages probably emerged or spread early in Portugal. We estimated a 62% infection reporting rate and that 1.3% of the population of men who have sex with men in Portugal were infected. We infer the critical role played by sexual networks and superspreader gatherings, such as sauna attendance, in the dissemination of mpox virus. Overall, our findings highlight genomic epidemiology as a tool for the real-time monitoring and control of mpox epidemics, and can guide future vaccine policy in a highly susceptible population.

X-Linked Myotubular Myopathy: A Novel Mutation Expanding the Genotypic Spectrum of a Phenotypically Heterogeneous Myopathy.

X-linked myotubular myopathy (XLMTM), a centronuclear congenital myopathy secondary to pathogenic variants in the MTM1 gene encoding myotubularin, is typically recognized for its classic and severe phenotype which includes neonatal hypotonia, severe muscle weakness, long-term ventilator dependence, markedly delayed gross motor milestones with inability to independently ambulate, and a high neonatal and childhood mortality. However, milder congenital forms of the condition and other phenotypes are recognized. We describe a 6-year-old boy with a mild XLMTM phenotype with independent gait and no respiratory insufficiency even in the neonatal period. The child has a hemizygous novel splice site variant in the MTM1 gene (c.232-25A > T) whose pathogenicity was confirmed by cDNA studies (exon 5 skipping) and muscle biopsy findings. We also compared the phenotype of our patient with the few reported cases that presented a mild XLMTM phenotype and no respiratory distress at birth, and discussed the potential mechanisms underlying this phenotype such as the presence of residual expression of the normal myotubularin transcript.

Muscle glycome in idiopathic inflammatory myopathies: Impact in IL-6 production and disease prognosis.

Idiopathic inflammatory myopathies (IIM) are a group of chronic autoimmune diseases mainly affecting proximal muscles. Absence of meaningful prognostic factors in IIM has hindered new therapies development. Glycans are essential molecules that regulate immunological tolerance and consequently the onset of autoreactive immune response. We showed that muscle biopsies from patients with IIM revealed a deficiency in the glycosylation pathway resulting in loss of branched N-glycans. At diagnosis, this glycosignature predicted disease relapse and treatment refractoriness. Peripheral CD4+ T cells from active-disease patients shown a deficiency in branched N-glycans, linked to increased IL-6 production. Glycan supplementation, restoring homeostatic glycosylation profile, led to a decrease in IL-6 levels. This study highlights the biological and clinical importance of glycosylation in IIM immunopathogenesis, providing a potential mechanism for IL-6 production. This pinpoints muscle glycome as promising biomarker for personalized follow-up and a potential target for new therapies in a patients' subgroup with an ominous evolution.

Looking at the Molecular Target of NS5A Inhibitors throughout a Population Highly Affected with Hepatitis C Virus.

Hepatitis C virus (HCV) is associated with liver damage and an increased progression rate to cirrhosis and hepatocellular carcinoma. In Portugal, it is prevalent in vulnerable populations such as injection drug users (IDU). HCV is characterized by a high intra-host variability, and the selecting driving forces could select variants containing resistance-associated substitutions (RAS) that reduce treatment effectiveness. The main goal of this study was to analyze the sequence variation of NS5A in treatment-naïve IDU. The epidemiological and clinical status of hepatitis C were analyzed, and samples were sequenced by Sanger and Next-Generation sequencing (NGS) to assess RAS and confirm HCV subtypes. Phylogenetic classification was concordant: 1a, 52.4%; 1b, 10.7%; 3a, 20.2%; 4a, 8.3%; 4d, 7.1%; and one 2k/1b recombinant. A 1a/3a mixed infection was detected by NGS. RAS were found in 34.5% (29/84) of samples using Sanger sequencing, while in 42.9% (36/84) using NGS. In sequences from subtypes 1a and 1b, RAS K24R, M28V, Q30H/R, H58D/P/Q/R, and RAS L31M and P58S were detected, respectively. In subtype 3a, RAS A30S/T, Y93H and polymorphisms in position 62 were identified. Additionally, RAS P58L was detected in genotype 4. The strategy used for the molecular survey of baseline HCV resistance is of particular importance to achieve treatment effectiveness and contribute to the elimination of hepatitis C.

ReporTree: a surveillance-oriented tool to strengthen the linkage between pathogen genetic clusters and epidemiological data.

BACKGROUND: Genomics-informed pathogen surveillance strengthens public health decision-making, playing an important role in infectious diseases' prevention and control. A pivotal outcome of genomics surveillance is the identification of pathogen genetic clusters and their characterization in terms of geotemporal spread or linkage to clinical and demographic data. This task often consists of the visual exploration of (large) phylogenetic trees and associated metadata, being time-consuming and difficult to reproduce. RESULTS: We developed ReporTree, a flexible bioinformatics pipeline that allows diving into the complexity of pathogen diversity to rapidly identify genetic clusters at any (or all) distance threshold(s) or cluster stability regions and to generate surveillance-oriented reports based on the available metadata, such as timespan, geography, or vaccination/clinical status. ReporTree is able to maintain cluster nomenclature in subsequent analyses and to generate a nomenclature code combining cluster information at different hierarchical levels, thus facilitating the active surveillance of clusters of interest. By handling several input formats and clustering methods, ReporTree is applicable to multiple pathogens, constituting a flexible resource that can be smoothly deployed in routine surveillance bioinformatics workflows with negligible computational and time costs. This is demonstrated through a comprehensive benchmarking of (i) the cg/wgMLST workflow with large datasets of four foodborne bacterial pathogens and (ii) the alignment-based SNP workflow with a large dataset of Mycobacterium tuberculosis. To further validate this tool, we reproduced a previous large-scale study on Neisseria gonorrhoeae, demonstrating how ReporTree is able to rapidly identify the main species genogroups and characterize them with key surveillance metadata, such as antibiotic resistance data. By providing examples for SARS-CoV-2 and the foodborne bacterial pathogen Listeria monocytogenes, we show how this tool is currently a useful asset in genomics-informed routine surveillance and outbreak detection of a wide variety of species. CONCLUSIONS: In summary, ReporTree is a pan-pathogen tool for automated and reproducible identification and characterization of genetic clusters that contributes to a sustainable and efficient public health genomics-informed pathogen surveillance. ReporTree is implemented in python 3.8 and is freely available at https://github.com/insapathogenomics/ReporTree .

Irreversible Geometrothermodynamics of Open Systems in Modified Gravity.

In this work, we explore the formalism of the irreversible thermodynamics of open systems and the possibility of gravitationally generated particle production in modified gravity. More specifically, we consider the scalar-tensor representation of f(R,T) gravity, in which the matter energy-momentum tensor is not conserved due to a nonminimal curvature-matter coupling. In the context of the irreversible thermodynamics of open systems, this non-conservation of the energy-momentum tensor can be interpreted as an irreversible flow of energy from the gravitational sector to the matter sector, which, in general, could result in particle creation. We obtain and discuss the expressions for the particle creation rate, the creation pressure, and the entropy and temperature evolutions. Applied together with the modified field equations of scalar-tensor f(R,T) gravity, the thermodynamics of open systems lead to a generalization of the ΛCDM cosmological paradigm, in which the particle creation rate and pressure are considered effectively as components of the cosmological fluid energy-momentum tensor. Thus, generally, modified theories of gravity in which these two quantities do not vanish provide a macroscopic phenomenological description of particle production in the cosmological fluid filling the Universe and also lead to the possibility of cosmological models that start from empty conditions and gradually build up matter and entropy.

Microevolution of a Mycobacteroides abscessus subsp. bolletii strain in a clinical persistent infection.

Mycobacteroides abscessus complex (MAB), a fast-growing nontuberculous mycobacterium, is emerging as a significant infectious disease threat, due to both intrinsic and acquired resistance mechanisms to antibiotics and disinfectants and the need for extensive and multidrug regimens for treatment. Despite the prolonged regimens, outcomes are poor and persistence cases have been reported. Here, we describe clinical, microbiologic and genomic features of a M. abscessus subsp. bolletii (M. bolletii) strain consecutively isolated from a patient within an eight-year infection period. From April 2014 to September 2021, the National Reference Laboratory for Mycobacteria received eight strains isolated from a male patient. Species identification, molecular resistance profile and phenotypic drug susceptibility were determined. Five of these isolates were recovered for further in-depth genomic analysis. Genomic analysis confirmed the multidrug resistant pattern of the strain and also other genetic changes associated with adaptation to environment and defence mechanisms. We highlight the identification of new mutations in locus MAB_1881c and in locus MAB_4099c (mps1 gene), already described as associated with macrolides resistance and morphotype switching, respectively. Additionally, we also observed the emergence and fixation of a mutation in locus MAB_0364c that appeared at a frequency of 36% for the 2014 isolate, 57% for the 2015 isolate and 100% for the 2017 and 2021 isolates, clearly illustrating a fixation process underlying a microevolution of the MAB strain within the patient. Altogether these results suggest that the observed genetic alterations are a reflection of the bacterial population's continuous adaptation and survival to the host environment during infection, contributing to persistence and treatment failure.

Comparative complete scheme and booster effectiveness of COVID-19 vaccines in preventing SARS-CoV-2 infections with SARS-CoV-2 Omicron (BA.1) and Delta (B.1.617.2) variants: A case-case study based on electronic health records.

Background: Information on vaccine effectiveness in a context of novel variants of concern (VOC) emergence is of key importance to inform public health policies. This study aimed to estimate a measure of comparative vaccine effectiveness between Omicron (BA.1) and Delta (B.1.617.2 and sub-lineages) VOC according to vaccination exposure (primary or booster). Methods: We developed a case-case study using data on RT-PCR SARS-CoV-2-positive cases notified in Portugal during Weeks 49-51, 2021. To obtain measure of comparative vaccine effectiveness, we compared the odds of vaccination in Omicron cases versus Delta using logistic regression adjusted for age group, sex, region, week of diagnosis, and laboratory of origin. Results: Higher odds of vaccination were observed in cases infected by Omicron VOC compared with Delta VOC cases for both complete primary vaccination (odds ratio [OR] = 2.1; 95% confidence interval [CI]: 1.8 to 2.4) and booster dose (OR = 5.2; 95% CI: 3.1 to 8.8), equivalent to reduction of vaccine effectiveness from 44.7% and 92.8%, observed against infection with Delta, to -6.0% (95% CI: 29.2% to 12.7%) and 62.7% (95% CI: 35.7% to 77.9%), observed against infection with Omicron, for complete primary vaccination and booster dose, respectively. Conclusion: Consistent reduction in vaccine-induced protection against infection with Omicron was observed. Complete primary vaccination may not be protective against SARS-CoV-2 infection in regions where Omicron variant is dominant.

Molecular investigation of ready-to-eat salads for Giardia duodenalis and Cryptosporidium spp. in Portugal.

Foodborne outbreaks are often associated with the consumption of salads. However, published studies on the detection of foodborne pathogens in ready-to-eat salads are scarce. The aim of this study was to detect Giardia duodenalis and Cryptosporidium DNA in ready-to-eat salads, by applying techniques of molecular biology to study the frequency of contamination in salads. A total of 100 packages of ready-to-eat salads containing assorted leafy green vegetables were randomly purchased from hypermarkets located in central regions of Portugal (Coimbra and Viseu). Nested-PCR and qPCR methods were used to detect G. duodenalis and Cryptosporidium DNA. Species and assemblages of the parasites were identified by sequence analysis and PCR. Eighteen of the 100 samples (18%) were positive for G. duodenalis and twelve were sequenced and identified as assemblage A. Cryptosporidium spp. were not detected in any salads. Overall, pre-harvest and post-harvest preventive measures may be need for G. duodenalis control throughout the food production industry, from the field to consumers.

Comparative Effectiveness of COVID-19 Vaccines in Preventing Infections and Disease Progression from SARS-CoV-2 Omicron BA.5 and BA.2, Portugal.

We estimated comparative primary and booster vaccine effectiveness (VE) of SARS-CoV-2 Omicron BA.5 and BA.2 lineages against infection and disease progression. During April-June 2022, we implemented a case-case and cohort study and classified lineages using whole-genome sequencing or spike gene target failure. For the case-case study, we estimated the adjusted odds ratios (aORs) of vaccination using a logistic regression. For the cohort study, we estimated VE against disease progression using a penalized logistic regression. We observed no reduced VE for primary (aOR 1.07 [95% CI 0.93-1.23]) or booster (aOR 0.96 [95% CI 0.84-1.09]) vaccination against BA.5 infection. Among BA.5 case-patients, booster VE against progression to hospitalization was lower than that among BA.2 case-patients (VE 77% [95% CI 49%-90%] vs. VE 93% [95% CI 86%-97%]). Although booster vaccination is less effective against BA.5 than against BA.2, it offers substantial protection against progression from BA.5 infection to severe disease.

Rapidly Progressive Corticobasal Degeneration Mimicking Brainstem Encephalitis.

Background: Corticobasal degeneration (CBD) may have a rapidly progressive (RP) clinical course, mimicking other neurological conditions. Objectives: To describe a neuropathologically proven case of RP-CBD in a patient initially diagnosed with immune-mediated brainstem encephalitis. Methods: Retrospective data collection from electronic records and authorized video material. Results: A 51-year-old man presented with bilateral ptosis, diplopia, and dysphagia. The diagnostic workup was negative for myasthenic syndromes. He progressively developed cognitive dysfunction with frontal release signs and asymmetric parkinsonism. Cerebrospinal fluid evaluation revealed 4 leukocytes/uL, 0.32 g/L proteins, 0.85 g/L glucose, and absent oligoclonal bands. Weakly positive anti-PNMA2 (Ma2/Ta) antibodies were present, and magnetic resonance imaging showed a T2 hyperintensity involving the midbrain and pons. Based on these features, the diagnosis of immune-mediated brainstem encephalitis was considered. The patient did not improve after several cycles of methylprednisolone, intravenous immunoglobulin, and plasma exchange. At 1 year after onset, he developed horizontal and vertical gaze limitation and worsening of the parkinsonism and cognitive dysfunction. By age 53, he was severely disabled, requiring percutaneous gastrostomy for feeding. Anti-IgLON5 was negative. He fulfilled the clinical criteria for probable progressive supranuclear palsy. He died from pneumonia at age 54. The neuropathological examination revealed a 4-repeat tauopathy with features of CBD with extensive involvement of the brainstem. Conclusions: RP-CBD may resemble brainstem encephalitis. The severity of brainstem and upper spinal cord pathology in the postmortem examination correlated with the clinical and imaging features.

Molecular Capture of Mycobacterium tuberculosis Genomes Directly from Clinical Samples: A Potential Backup Approach for Epidemiological and Drug Susceptibility Inferences.

The application of whole genome sequencing of Mycobacterium tuberculosis directly on clinical samples has been investigated as a means to avoid the time-consuming need for culture isolation that can lead to a potential prolonged suboptimal antibiotic treatment. We aimed to provide a proof-of-concept regarding the application of the molecular capture of M. tuberculosis genomes directly from positive sputum samples as an approach for epidemiological and drug susceptibility predictions. Smear-positive sputum samples (n = 100) were subjected to the SureSelectXT HS Target Enrichment protocol (Agilent Technologies, Santa Clara, CA, USA) and whole-genome sequencing analysis. A higher number of reads on target were obtained for higher smear grades samples (i.e., 3+ followed by 2+). Moreover, 37 out of 100 samples showed ≥90% of the reference genome covered with at least 10-fold depth of coverage (27, 9, and 1 samples were 3+, 2+, and 1+, respectively). Regarding drug-resistance/susceptibility prediction, for 42 samples, ≥90% of the >9000 hits that are surveyed by TB-profiler were detected. Our results demonstrated that M. tuberculosis genome capture and sequencing directly from clinical samples constitute a potential valid backup approach for phylogenetic inferences and resistance prediction, essentially in settings when culture is not routinely performed or for samples that fail to grow.

Neuropathology of central nervous system involvement in TTR amyloidosis.

Hereditary transthyretin amyloidosis (ATTRv) is a systemic disease caused by the accumulation of misfolded transthyretin (TTR). It usually presents with an adult-onset progressive axonal peripheral neuropathy and cardiomyopathy. In the central nervous system (CNS), variant TTR is produced by the choroid plexus and accumulates in the leptomeninges. CNS symptoms have been increasingly recognized in this population, including transient focal neurological episodes and stroke, particularly in patients with the V30M mutation and longstanding disease. The prevalence, pathophysiology, and progression of CNS involvement remain to be clarified. The present work explores if there is a recognizable sequence of CNS TTR deposition in ATTRv. We studied the topographical and severity distribution of TTR deposition in 16 patients with ATTRv, aged 27-69 years and with a mean disease duration of 10.9 years (range: 3-29). Our results suggest that CNS pathological involvement in V30M ATTRv occurs early in the disease course, probably starting in pre-symptomatic phases, and follows a distinct sequence. Leptomeninges and subarachnoid meningeal vessels are affected earlier, then followed by perforating cortical vessels and subpial deposition, and finally by deposition in the subependymal and basal ganglia vessels near the ependymal lining. Brainstem and spinal cord show early and severe involvement, with amyloid subpial deposition already seen in initial stages. Despite massive superficial amyloid deposition, no parenchymal deposition outside subpial or subependymal regions was found. Additionally, vascular lesions or superficial cortical siderosis were not frequent. Future studies with more patients from different populations and TTR mutations will be important to confirm these findings. Defining stages of TTR pathology in the CNS may be useful to better understand pathogenic mechanisms leading to symptoms and to interpret neuroimaging biomarkers.

Strengthening the genomic surveillance of Francisella tularensis by using culture-free whole-genome sequencing from biological samples.

Introduction: Francisella tularensis is a highly infectious bacterium that causes the zoonotic disease tularemia. The development of genotyping methods, especially those based on whole-genome sequencing (WGS), has recently increased the knowledge on the epidemiology of this disease. However, due to the difficulties associated with the growth and isolation of this fastidious pathogen in culture, the availability of strains and subsequently WGS data is still limited. Methods: To surpass these constraints, we aimed to implement a culture-free approach to capture and sequence F. tularensis genomes directly from complex samples. Biological samples obtained from 50 common voles and 13 Iberian hares collected in Spain were confirmed as positive for F. tularensis subsp. holarctica and subjected to a WGS target capture and enrichment protocol, using RNA oligonucleotide baits designed to cover F. tularensis genomic diversity. Results: We obtained full genome sequences of F. tularensis from 13 animals (20.6%), two of which had mixed infections with distinct genotypes, and achieved a higher success rate when compared with culture-dependent WGS (only successful for two animals). The new genomes belonged to different clades commonly identified in Europe (B.49, B.51 and B.262) and subclades. Despite being phylogenetically closely related to other genomes from Spain, the detected clusters were often found in other countries. A comprehensive phylogenetic analysis, integrating 599 F. tularensis subsp. holarctica genomes, showed that most (sub)clades are found in both humans and animals and that closely related strains are found in different, and often geographically distant, countries. Discussion: Overall, we show that the implemented culture-free WGS methodology yields timely, complete and high-quality genomic data of F. tularensis, being a highly valuable approach to promote and potentiate the genomic surveillance of F. tularensis and ultimately increase the knowledge on the genomics, ecology and epidemiology of this highly infectious pathogen.